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Abstract 

' Radiative decays of the l^Si radial and l^Di orbital excitations of the p, uj and (j> are calculated 

in the quark model, using wave functions obtained variationally from the Hamiltonian with 
standard quark-model parameters. The larger radiative widths should be measurable at new 
high-intensity facilities being proposed, and in some cases may be measurable in data from 
existing experiments. The radiative decays are a strong discriminator between the l^Si and 
^ • l^Di excitations, and can also be used to provide unique information about the decay products. 
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1 Introduction 



^ , Radiative decays offer a rather direct probe of hadron structure. The coupling to the charges 

' and spins of constituents reveals detailed information about wave functions and can discriminate 

among models. This can be particularly relevant in distinguishing gluonic excitations of the tt, p 
(hybrids) from conventional excitations with these overall J^^ . For example, in a hybrid 1 the 
qq are in a spin-singlet, while for the ^ they are in a triplet: in each case this is the reverse of 
what one is used to. 

In this paper we investigate potential tests using radiative decays of vector mesons to separate 
qq states from hybrids in the 1-2 GeV mass region, where light-flavour states are predicted to 
occur, and where prima facie candidates have been identified. Exotic J^^ — 1 ^ signals have 
been reported around 1.4-1.6 GeV mass, the systematics of 1 states in this region seem to point 
towards gluonic excitations being present and there are tantalising hints of unusual activity in the 
^ partial wave in the 1.6-1.8 GeV region. 

A quite separate problem is posed by the scalar mesons. It seems likely that the ao(980) and 
/o(980) are intimately linked to the KK threshold with significant q^if affinity. Then the qq 
mesons are manifested in the 1.3 to 1.7 GeV region, presumably mixed with the 0^'*' glueball of 
the lattice. A direct measure of their flavour content would resolve this issue. We suggest that 
radiative transitions from the 25* and ID vector excitations can answer this. 

The forthcoming generation of e+e^ facilities, such as the upgrade of VEPP at Novosibirsk, the use 
of initial state radiation (ISR) at BABAR to study light-quark vectors and the proposed dedicated 
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e+e~ collider PEP-N, promise more than two orders of magnitude increase in data over present 
machines. There is a complementary high-intensity programme of photo- and electroproduction at 
Jefferson Laboratory. Together these will enable data on radiative decays to be obtained compa- 
rable in quality to that of present hadronic channels. Thus it is timely to ask whether, and under 
what circumstances, these data will be able to make sensitive probes of radiative transitions and to 
isolate evidence for gluonic degrees of freedom. Radiative decays offer an alternative approach to 
resolving the various problems we have outlined. They are a much better probe of wave functions, 
and hence of models, than are hadronic decays. They can access final states which are kinemat- 
ically excluded for hadronic decays. Historically, studies of light-quark radiative transitions have 
been restricted to the ground states. We extend these calculations to decays from the 25" and 
ID excitations of the p, uj and to 15, 25' and IP states, as these processes are now becoming 
accessible to experiment. 

Our results are encouraging. We find that certain channels, with large partial widths, allow a clean 
separation of the 25 and \D qq states, for example /2(1270)7 for the /3(1450), /17 for the p(1700), 
/;^(1525)7 for the (/i(1690) and /i(1420) for the (/)(1900). Additionally the radiative decays of the 
/9(1700) and (/)(1900) to nn and ss scalars respectively, probe the /o sector and provide a potential 
entree to flavour filtering. 

In section 2 we review the present status of the interesting mesons in the 1 to 2 GeV mass range. 
We then turn in sections 3 and 4 to establishing our formalism for calculating radiative transition 
rates. Their magnitudes are computed in section 5, and in section 6 we identify a strategy for 
exploiting these results in forthcoming experiments. 

2 Mesons in the 1 to 2 GeV mass range 

The existence of two higher isovector vector mesons, the p(1450) and the p(1700); their isoscalar 
counterparts, the cl'(1420) and a;(1650); and of an associated hidden-strangeness state, the (^(1680), 
is well established Although there is general consensus on the existence of these states, there is 
considerable disparity on their masses and widths. Further, what is known about the composition 
of their hadronic decays raises fundamental questions about the nature of these states and our 
understanding of the mechanism of hadronic decays. 

An apparently natural explanation for the higher-mass vector states is that they are the first radial, 
2'^5i, and first orbital, l^Di, excitations of the p and uj and the first radial excitation of the cf), 
as the generally-accepted masses jl) are close to those predicted by the quark model [1). However 
this argument is suspect as the masses of the corresponding ~ 1^ strange mesons are less than 
the predictions, particularly for the 2'^5i at 1414 ± 15 MeV [|| compared to the predicted 1580 
MeV 1^. Quite apart from comparing predicted and observed masses, one would expect the nfi 
mesons to be 100 to 150 MeV lighter than their strange counterparts, putting the 2'^5i at less 
than 1300 MeV and the I'^Di below 1600 MeV. Also this interpretation faces a further problem. 
The data on the 47r channels in e~^e~ annihilation and r decay are not compatible with the ^Pq 
model |, § |, |, which works well for decays of established ground-state mesons. For example 
widths predicted to be large are found to be so; widths predicted to be small, are found to be so; 
calculated widths agree with data to 25 — 40%; and signs of amplitudes are predicted correctly. 
As far as one can ascertain the '^Pq model is reliable, but it has not been seriously tested for the 
decays of excited states. 

The '^Pq model predicts that the decay of the isovector 2'^5i to An is extremely small: r2s^ai7r ~ 3 
MeV and r2s^/ii7r ~ 1 MeV, and other possible An decays are even smaller. For the isovector I'^Di 
the model predicts that the ain and hin decays are large and equal: Tiu^aiTr ^ ^iD^hiir ^ 105 
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MeV. All other possible 47r channels are small. As hiir contributes only to the 7r"'"7r~7r''7r° channel in 
e~^e~ annihilation, and oitt contributes to both 7r+7r~7r"'"7r~ and 7r"'"7r~7r°7r'', then after subtraction 
of the u!Tr cross section from the total 7r+7r^7r'^7r'^, one expects that cr(e+e^ — > 7r+7r~7r°7r°) > 
cr(e"'"e~ Tr+TT^TT+TT^). This contradicts observation 0| over most of the available energy range, 
in which cr(7r+7r^7r+7r^) « 2cr(7r"'"7r~7r'^7r°). Further, and more seriously, it has been shown recently 
by the CMD collaboration at Novosibirsk 0] and by CLEO that the dominant channel by far 
in Att (excluding lott) up to 1.6 GeV is aiir. This is quite inexplicable in terms of the ^Pq model. 
So the standard picture is wrong for the isovectors, and there are serious inconsistencies in the 
isoscalar channels as well. One possibility is that the ^Pq model is simply failing when applied to 
excited states, which is an intriguing question in itself. An alternative is that there is new physics 
involved. 



A favoured hypothesis is to include vector hybrids ||9|, that is qqg states. The reason for this 

is that, firstly, hybrid states occur naturally in QCD, and secondly, that in the relevant mass range 
the dominant hadronic decay of the isovector vector hybrid pn is believed to be oitt [|io||. The 



masses of light-quark hybrids have been obtained in lattice-QCD calculations ||T^, |l^, 14, [15|, 
although with quite large errors. Results from lattice QCD and other approaches, such as the bag 
model [|l^ |l^, flux-tube models ||l9|, constituent gluon models |^ and QCD sum rules [|T| p^ , 
show considerable variation from each other. So the absolute mass scale is somewhat imprecise, 
predictions for the lightest hybrid lying between 1.3 and 1.9 GeV. However most models and lattice 
gauge calculations predict the lightest hybrids to be at the upper end of this mass range. It does 
seem generally agreed that the mass ordering is < 1 < 1 < 2 ^. 

Evidence for the excitation of gluonic degrees of freedom has emerged in several processes. A clear 
exotic J^'-^ = 1 ^ resonance, the 7ri(1600), is seen in the 77'(958)7r channel in the reaction 
TT" N — > (?7'(958))7riV. Two experiments |2^, have evidence for this exotic in the p'^Tr" channel 
in the reaction ir^N — > (7r+7r^7r~)A^. A peak in the rjir mass spectrum at ^ 1400 MeV with 
JP^ - 1- in TT N ^ (ryTT )N has also been interpreted as a resonance [^. Supporting evidence 
for the 1400 MeV state in the same mode comes from pp rjir'TT^ There is evidence for 
two isovector ^ states in the mass region 1.4 to 1.9 GeV; 7r(1600) and 7r(1800). The quark model 
predicts only one. Taking the mass of the 1 ^ 1.4 GeV and assuming the generally-agreed mass 
ordering, then the is at ~ 1.3 GeV and the lightest 1^~ at ~ 1.65 GeV, which is in the range 
required for the mixing hypothesis to work. However if the 7ri(1600) is the lightest 1 ^ state, 
so that the vector hybrids are comparatively heavy, say ~ 2.0 GeV, then strong mixing with the 
radial and orbital excitations would be less likely. At this stage we should keep an open mind and 
consider both options. 

The scalar mesons in the mass range 1.3 to 1.7 GeV provide another place in which to look for 
gluonic hadrons as lattice calculations have now become sufficiently stable to predict |2|], in the 
quenched approximation, that the lightest glueball has J^'-^ = 0"'"'" and is in the mass range 1.45 
to 1.75 GeV. 

There are more scalar mesons than the simple qq l^Po nonet can accomodate. The /o(980) and 
ao(980) mesons most probably do not belong to this nonet |3^, and are either q^(f states or 
KK molecules. In any case they are associated with the nearby KK threshold. Then the possible 
candidates for the qq l^Pg nonet are ao(1450), K^{Um), /o(1370), /o(1500) and /o(1710). There is 
an obvious excess in the isoscalar-scalar sector, with the natural inference of there being a glueball 
state present. 

The lightest scalar glueball should mix with the qq scalars in the same mass region and recent 
studies J32] on a coarse-grained lattice suggest that such mixing is significant. While analyses 
pO| , |3^ , |34|, |35[| of the mixing differ in some details, the conclusions exhibit common robust features. 
The flavour content is predicted to have nfi and ss in phase for the /o(1370) and /o(1710) (SU(3)- 
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singlet tendency), out of phase for the /o(1500) (SU(3)-octet tendency) and to have a glueball 
component in all three states. The detailed pattern of mixing was determined in Jssf by studying 
the complete set of decay branching ratios into pseudoscalars Q for the /o(1370), /o(1500) and 
/o(1710), and confirmed by comparing relative production rates. The preferred scenario Q gives 
the bare masses as rUg = 1443 ± 24 MeV, mnn = 1377 ± 10 MeV and mss = 1674 ± 10 MeV. 
Other solutions have been found which have either a heavy glueball, rUg > russ, or a light glueball, 
■nig < rrinn, and although less consistent with the data they cannot be ruled out completely. The 
preferred solution is consistent with what one would expect naively from the ss — nn mass difference 
of about 300 MeV, and places the glueball at the lower end of the mass range given by the lattice 
calculations. 

The mixing scheme implies that the isovector partner of the nn state, the uq, should have a mass of 
about 1400 MeV. There is an indication that this state has been observed ||3^. Any confirmation 
of this controversial ao(1450) is of paramount importance, not only for the problem of glueball- 
quarkonia mixing, but for the nature of the 0^^ mesons. 

The emerging data suggest that gluonic excitations, both for hybrids and glueballs, are rather 
lighter than quenched-lattice predictions and that their effects will be apparent in the 1 to 2 GeV 
mass range. As we shall see, radiative transitions can shed new light on the matter. 



3 Radiative Decays of Quarkonia 

The initial meson A, with mass ttia, decays at rest to the final-state meson B, with mass and 
a photon with three-momentum p. In the non-relativistic quark model the standard expression for 
the transition amplitude has the form 

Ma^b = M^^5 + M\^^ (1) 

where M^_^^ and M^^^ describe the emission of the photon from the quark and antiquark 
respectively: 

-irr{</.^3(k-ip)a0^(k)}xp] (2) 

and 

Mi_B - /d3fc[Tr{0^(k)0t^(k+ip)}(2k + p) 

Zniq J 

-tTr{M^)acpl{k +^p)}xp] (3) 
where Iq and Iq are isospin factors and mq is the quark mass. 

We use matrix forms for the wave functions. For a meson M, with quark spin 0, total angular 
momentum j and magnetic quantum number m the wave function is given by 

1 



0M(q) = —iY,^{q)RM{q) ■ (4) 



For a meson with total angular momentum j, quark spin 1 and quark orbital momentum /, the 
corresponding wave function is 

</'Af(q) = ■^=Y jim{q)cTRM{q) ■ (5) 
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Here 1 is the 2x2 unit matrix and a is the PauU matrix. The Ruifl) arc the mesonic radial wave 
functions in the momentum representation 

We calculate in the centre-of-mass of e^e~ annihilation as in that case the virtual photon is 
polarised, and one can study angular distributions with respect to the beam direction. These are 
given in Appendix A. In Appendix B we show how to derive the helicity amplitudes and their 
relation to the formalism used in the paper, which enables the results to be transferred to any 
frame of reference. 

The differential decay rate is evaluated for initial photon polarisation toq = 1 and is given by 

-^^^^V—o^lY.\^lA^^? (6) 
dcoso niA ^ 

where the sum is over final-state polarisations. In Eb is the centre-of-mass energy of the final 
meson and I = Ig = is the isospin factor. We consider the radiative decays of neutral vector 
mesons, so the isospin factors for decays between nfi or ss states are 

/ = for nn nn with same isospin 

I = J for nn nn with different isospin 

/ = ^ for ss —> ss. (7) 

We take rj — rjs ~ {uu + dd— 2ss) / and rj' = rji — {uu + dd + ss) / \/3 so the isospin factors for 
decays to rj and rj' are 

/ = jig for nn with isospin — > 77 

I = j2 for with isospin 1 ^ 77 

/ = for nn with isospin ~-^ rj' 

/ = i for nn with isospin 1 ^ rj' 

I = -§f for ss — > ?y 

/ = ^ for ss 77'. (8) 



Radial wave functions are found variationally from the Hamiltonian 

H=^+ar-^^ + C (9) 
niq 6 r 

with standard quark model parameters mq = 0.33 GeV for u and d quarks, niq = 0.45 GeV for s 
quarks, a = 0.18 GeV^ and as = 0.5. The wave functions are taken to be Gaussian, that is of the 
form exp(— A:^/(2/3^^)) multiplied by the appropriate polynomials, and /3 treated as the variational 
parameter in (^) for each of the IS", IP, 2S and \D states. The resulting values of (3 and the 
corresponding masses for nn states are given in table 1, and those for ss states in table 2. 



It is well-known that the use of exact wave functions is necessary to reproduce the low-energy 
theorems. For example, the expression for the El transition amplitude (equations and in 
the p — Q limit) can be written in the dipole form using the relation k = innq[Hr\/2. The amplitude 
is proportional to the overlap integral 

'■^{Ea - Eb) I r^drijA{r)iJB{r), (10) 
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M{GeV) 


f3{GeV) 


IS 


0.700 


0.313 


IP 


1.262 


0.274 


2S 


1.563 


0.253 


ID 


1.703 


0.255 



Table 1. Effective masses and corresponding /? from the variational solution of (|9|) for nn states. 





M{GeV) 


f3{GeV) 


IS 


1.000 


0.355 


IP 


1.527 


0.307 


2S 


1.793 


0.285 


ID 


1.932 


0.285 



Table 2. Effective masses and corresponding f3 from the variational solution of (|9|) for ss states. 



which immediately gives the threshold behaviour of the electromagnetic width demanded by 
gauge invariance. This means that not only should the "correct" values of the oscillator parameter 
(3 be used, but also the masses of the initial and final states should be taken to be "correct" 
eigenvalues of the quark model Hamiltonian (^). One can systematically improve the variational 
ansatz adopted for the wave functions, using more sophisticated trial wave functions ||^, |3^ or 
numerical solutions |Q. Fortunately, in most cases in practice, the "correct" masses and the 
known physical masses do not differ by much, and the resulting differences in the decay widths 
are smaller than other uncertainties in the calculation. However this is not true of decays to the 
pseudoscalars tt, rj and r/'. 

The constituent quark model in its naive form works reasonably well, with the exception of the 
lowest pseudoscalars as the Goldstone-boson nature of these particles is not naturally accomodated. 
More complicated versions, such as ||^, include relativistic kinematics and properly-smeared spin- 
dependent terms, and the agreement with data on spectra is quite remarkable. Similar results are 
obtained in the QCD string model ||3^, where it is shown that the constituent masses of quarks 
appear dynamically due to the QCD-inspired string-type interaction. The most serious drawback 
of such a picture is that the issue of chiral symmetry breaking is completely beyond its scope. 
It is possible to describe the low-lying pion in the constituent quark models, but the Goldstone 
nature of the pion is completely lost. However, the chiral properties can be naturally incorporated 
into the constituent picture in the framework of the Hamiltonian approach in the Coulomb gauge 
. In this approach, a chirally-noninvariant vacuum is constructed, which implies the existence 
of the Goldstone boson as the lowest qq state. The axial current is conserved in the chiral limit, 
and all the relations of current algebra are satisfied. Some results of this approach are relevant 
to our purposes. First, even in the chiral limit, where the pionic wave function has very peculiar 
properties, the other mesons behave to large extent as quark model ones (for the details see 
Second, while the confining interaction alone can describe the observed n-p splitting, the predictions 
for the quark condensate and pion decay constant are too small, and can be improved by inclusion 
of the hyperfine interaction Q. The net result of these studies is that if the effective degrees 
of freedom are properly defined, then the constituent quark model is a good approximation, with 
most important ingredients absorbed in the quark model parameters. 

This means that, as a first approximation, the mesonic wave functions can be found from the 
Hamiltonian (^. Also, for the pion, the "correct" (3 should be larger than the 313 MeV for 15 
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states shown in table 1. We will return to this point below. 



4 Radiative Decay Widths of Vector Mesons 

Expressions for the full angular distributions for the radiative decay widths of the I'^S'i, and 
l^Di neutral vector mesons are given in Appendix A. Here we give only the total widths. In these 
results we define 

^ {Pl + f3l) ^ ' 

and ^ 

A method for deriving the following expressions when /3a = Pb is given in Appendix B. The 
generalisation to arbitrary (3a and f3B is straightforward but algebraically tedious. 

l^Si l^So 

Q 4 Eb 2 ni\ 



2^Si l^So 



2^81 ^ 2^So 



4VV3/3l/3^ J 3/32 \3f3j,(3i, I 



(3M^"^) + 9mIi ^^^^ 



In (13), mj) and (hj Fsi is defined by 



23Si ^ l3Po 



23Si ^ l^Pi 



In (|17|), (|l|) and @ is defined by 
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Ge by 
and Gm by 

l^Di ^ l^So 

l^Di ^ 2^So 



G 



^ + (^-l)(6A-3).^(2A-l). (22) 



In (|23|) and (gj) is defined by 



l^Di ^ l^Pj 



^^^rA3(-2 + 3A + 2A2) + ^X\l - 4A + 4A2))/f|,2, (27) 



27 niA m-t V 5p^ 



+ YU^A'(83 - 316A + 520A2) + - 6A + 8A^) 



+ ^A^(1-4A + 4A^))/F^. (28) 



In (26), (M) and (E8h F^a is defined by 



5 Results 



We first give the numerical results for all relevant radiative decays. Discussion of these results is 
in two parts. The first treats decays to ^Pj states, with the emphasis on unique signatures for 
specific vector qq states. The second deals with radiative decays to scalars as a flavour filter with 
implications for the scalar glueball mass. We then discuss qualitatively hybrid meson radiative 
transitions. 
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r, 



700 



Lj — >7r7 

(keV) 



600 




300 - 
200 - 
100 - 



I ' ' ' ' ' 1 

0.2 0.25 0.3 0.35 0.4 0.45 0.5 

P. (GeV) 



Figure 1: Variation of the width for the decay ^ 7r7 as a function of with (3^^ fixed at 0.313 
according to table 1. 



Decay 


Model 


Experiment 




39.6 


36± 13 


w — > 777 


4.7 


5.5 ±0.9 




92.7 


57.82 ± 1.53 


(j) 


0.35 


0.30 ±0.15 



Table 3. Radiative widths for l^Si decays to 1^5*0. 



5.1 Numerical results 

A simple test of the vahdity of tables 1 and 2 is provided by the well-known decays p rjj, 
u> 777, (/) ^ ?77 and rj'j. The calculated widths in keV, using /3 = 0.313 from table 1 for 
the nn decays and /? = 0.355 from table 2 for the ss decays but the physical masses in each case, 
are compared with the experimental widths ^ in table 3. The agreement for the p and lu decays 
is clearly satisfactory, that for the decay (j) ^ -q^ less so. 

We can use the decays p — s- tt^ and lj — > 7r7, again with physical masses, to estimate the appropriate 
value of /Jtt for decays to tt^. The variation of the width for w — > 7r7 with (3t, is shown in figure 1. 
At no point does it reach the experimental value [|] of 717 ± 43 keV. The maximum of the curve 
is 520 keV and occurs at — 0.335, which we take to be the optimum value. The discrepancy 
between the model calculation and the experimental value is an indication of the uncertainty in 
the evaluation of decays to t:^. 

The radiative decay widths for the nn states are given in table 4. For purposes of illustration, we 
have assumed that the /o(1370) is a pure nn state and have ignored possible mixing in the I'^Pq 
nonet. This is certainly an oversimplification as there is good evidence discussed in section 2, 
for scalar mixing within this nonet and with a scalar glueball. We return to this question below. 
The effect of phase space is clearly seen in the deviations from the naive 9:1 ratios of (|^) and (||). 
The radiative decay widths for the ss states are given in table 5. As for the nn decays, we have 
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r(p(1450)) 


r(w(1420)) 


r(p(1700)) 


r(w(1650)) 




61 


510 


3.8 


40 


?77 


106 


11 


12 


1.1 


?7'7 


61 


5.7 


6 


0.5 


7r(1300)7 


5.9 


29 


0.4 


1.7 


77(1295)7 


57 


3.6 


3.9 


0.2 


/o (1370)7 


64 


4.7 


899 


88 


ao(1450)7 






82 


612 


/i(1285)7 


349 


33 


1097 


106 


ai(1260)7 


43 


341 


129 


1016 


/2 (1270)7 


712 


67 


148 


13 


02(1320)7 


59 


413 


13 


91 



Table 4. Results for radiative decays of nn states in keV. 





r((/)(1680)) 


r((/)(1930)) 


777 


94 


9 


v'l 


21 


1.8 


77(1440)7 


47 


10 


./o (1710)7 




188 


/i (1420)7 


148 


408 


/2(1525)7 


199 


37 



Table 5. Results for radiative decays of ss states in keV. 



ignored mixing and assumed that the /o(1710) is pure ss. 

We have commented above on the uncertainties in the decays p ^ tt'j and w — > 7r7. The situation 
becomes even more uncertain for the decays p{lA50) and w(1420) — > 7r7. These decays 

proceed essentially via the spin-flip part of the amplitude which vanishes in the nonrelativistic 
limit for orthogonal wave functions. For Gaussian wave functions the amplitude is proportional to 
{(3'^ / (3\ — 1) + 2A^p^/(3/9^), see (|4|), where the term {(3'^ / (3\ — 1) measures the nonorthogonality 
of the wave functions. Note that one should not expect the radial wave functions of the 2^iS'i 
and 1^5*0 states to be orthogonal in the fully relativistic theory, as spin-dependent forces are not 
treated as perturbations there. With f^A ^ Pb the resulting width is rather sensitive to the actual 
value of (3b ■ The value of /Jtt is not determined by the model, and as the width for lo tt^ varies 
rather slowly with Pt^, our choice of = 0.335 has a rather large error. The effect on the width 
of w(1420) — !■ 7r7 of varying is shown in figure ^. 

The first term in (|lj) increases rapidly with increasing /J^r , and for in the range we are considering 
it is this term and the interference term which dominate the width. The direct contribution from 
the second term in (|lj) becomes increasingly unimportant. A similar result holds for the decays 
p(1450) rjj and aj(1429) — > 777 as can be seen in figure ||. As a consequence we suggest that the 
errors on the 25* decays to ttj and r;7 are of the order of ±50%. 

An estimate of the uncertainties for other decays can be obtained by looking at the contributions 
from terms in increasing powers, n, of p^//3^. These are given in table 6, as percentages, for some 
of the larger widths. 
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800 



- iij(1420)^7r7 



(keV) 



600 



400 - 



200 - 



-200 



0.24 




0.32 0.34 0.36 

P. (GeV) 



Figure 2: Variation of the width for the decay ti>(1420) ^ 777 as a function of (3.^. The sohd hne is 
the total width, the dashed hne is the contribution from the first term in ([T^), the dotted hne is 
the contribution from the second term in (|lj) and the short-dashed hne is the contribution from 
the interference term. 



200 



p(1450)->r;7 



(keV) 



150 



100 - 




0.36 



(GeV) 



Figure 3: Variation of the width for the decay p(1450) ^ 777 as a function of The sohd hne is 
the total width, the dashed line is the contribution from the first term in ([T^), the dotted line is 
the contribution from the second term in ( |l^ ) and the short-dashed line is the contribution from 
the interference term. 
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n=0 


n=l 


n=2 


n=3 


n=4 


p(1450) ^ /i(1285)7 


106.2 


-6.4 


0.2 






p(1450) ^ /2(1270)7 


93.1 


6.6 


0.3 






p(1700) ^ /o(1370)7 


85.4 


14.8 


-0.2 


-0.1 


0.1 


p(1700) ^ /i(1285)7 


56.5 


39.6 


4.7 


-0.9 


0.1 


(/)(1680) ^ /i (1420)7 


112.7 


-13.5 


0.8 






(/)(1680) ^ /2(1525)7 


95.0 


4.8 


0.2 






(/.(1900) ^ /o(1710)7 


94.7 


5.3 


0.0 


0.0 


0.0 


0(1900) ^ /i(1420)7 


52.1 


43.2 


6.0 


-1.4 


0.1 



Table 6. Contributions, in percent, from terms in increasing powers, n, of p^//3^. 



The implications of table 6 are that the results for the decays of /5(1450) to /i (1285)7 /2( 1270)7, 
and the corresponding a;(1420) decays to ai(1260)7 and 02(1320)7 are very reliable, the decay of 
p(1700) to /o(1370)7 is solid, but the decay of p(1700) to /i(1285)7, and the corresponding w(1650) 
decay to ai (1260)7 should be treated with caution. This is not suprising because of the considerable 
increase in phase space. The situation is similar for the decays of 0(1690) and 0(1900), with the 
/i (1420)7 decay of the latter being the only one with any degree of uncertainty. 



5.2 Decays to ^Fj states 

It is clear from table 4 that some decays provide very clear signatures for particular excitations. 
Obvious ones are p{M50) /27 and aj(1420) 027; /9(1700) /07 and cj(1650) — > 097; 
/9(1700) —^ /i7 and tj(1650) —^ 017. However, in e"*"e~ annihilation experiments, isovector states 
are produced at approximately nine times the rate of the corresponding isoscalar states. Thus, for 
example, the effective 007 rate from the decay of tj(1650) will be the same as that from p(1700). 
Conversely the contamination of the /27 decay of /9(1450) by the corresponding decay of a;(1420) 
will only be at the 1% level. The /2 (1525)7 decay of the 0(1690) is important as it provides a 
unique signature for the ss state, since in contrast to hadronic decays it is unaffected by w(1650) 
and 0(1690) mixing. As we argue in section 5.4 radiative decays of a hybrid pn to /27 and of a 
hybrid (pn to /2 (1525)7 are strongly suppressed, so there is no ambiguity. 

Although the /i (1420)7 decay of the 0(1900) looks like a clean signature of this state, it should 
be recalled that the value quoted in table 5 should be treated with some caution. Despite this 
uncertainty, the width is necessarily much larger than that for the ry( 1440)7 decay. This provides 
a mechanism for producing the /i(1420) without contamination from 77(1440) with which it shares 
many common hadronic decay channels. 

It should be noted that there is some uncertainty about the mass of the p(1450) and aj(1420). For 
the isovector states, the most extreme low mass comes from an analysis of of the tt+tt^ spectrum 



in the reaction K~p 7r+7r~A |43 , which gives 1266 ± 14 MeV. An equally low mass, 1250 ± 29 
MeV has been suggested [Q for the isoscalar channel from an analysis of e+e~ tt+tt^tt^. Given 
this uncertainty, we show the mass variation of the width for the decay of the isovector radial 
(2'^ 5*1) excitation to /27 in figure 2. 



5.3 Scalar mesons and glueballs 

In tables 4 and 5 we assumed that there is no mixing among the scalars, so that the /o(1370) is 
pure nn and the /o(1710) is pure ss. The result of the mixing is that the bare nn and ss states 
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Figure 4: Variation of the width for the decay p{2S) — > /27 as a function of the mass of the p{2S). 

contribute in varying degrees to each of the /o(1370), /o(1500), and /o(1710). The variation of the 
radiative decay width of the p(1700) and 0(1900) as a function of the mass of the /o are shown in 
figures ^ and || respectively. 

Three different mixing scenarios have been proposed: the bare gluebaU is hghter than the bare nn 
state (the hght gluebaU solution); the mass of the bare glueball is between the bare nn state and 
the bare ss state (the middleweight glueball solution); and the mass of the bare glueball is greater 
than the mass of the bare ss state. The first two solutions have been obtained in |35| and the 
third has been suggested in iQ . The effects of the mixing on the radiative decay widths of the 
/9(1700) and the (/)(1900) to the three /q states are given in table 7 for each of these three cases. 
The relative rates of the radiative decays of the /9(1700) to /o(1370) and /o(1500) change radically 
according to the presence of the glueball admixture. So for a light glueball the decay to /o(1370) 
is relatively suppressed whereas for a heavy glueball it is substantial. By contrast the effect on 
the decay to /o(1500) goes the other way. Further, the (/)(1900) would give a large width for the 
decay to /o(1500) for a heavy glueball, but essentially zero for a light one. The /o(1710) will be 
prominent in the decays of the 0(1900) for all but the heaviest glueball. It is clear that these 
decays do provide an effective flavour-filtering mechanism. 





L 


p(1700) 
M 


H 


L 


0(1900) 
M 


H 


/o(1370) 


174 


440 


603 


7 


8 


31 


/o(1500) 


520 


301 


98 


5 


35 


261 


/o(1710) 




173 


156 


17 



Table 7. Effect of mixing in the scalar sector of the I^Pq nonet. The radiative widths, in keV, are 
given for three different mixing scenarios as described in the text: light glueball (L), middleweight 
glueball (M) and heavy glueball (H). 

Further, identifying the appropriate mixing scheme gives insight into the underlying physics of 
glueballs. The existing phenomenology from hadronic decays seems to favour a light glueball. 
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Figure 5: Variation of the width for the decay /9(1700) /07 as a function of the mass of the /q. 
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Figure 6: Variation of the width for the decay 0(1900) /07 as a function of the mass of the /q. 
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Essentially, if the decays of the "bare" glueball are flavour-independent, then the observed flavour 
dependencies for the hadronic decays of the physical mesons require jS^ the glueball mass to be 
at the low end of the range preferred by qucnchcd-lattice studies. 

The resolution of the isoscalar-scalar problem is intimately connected with the isovector-scalar 
problem. The existence of any aq other than the ao(980) remains controversial. The different 
mixing schemes for the isoscalar-scalar mesons give rather different values for the mass of the bare 
nfi state. This mass will be reflected in the mass of its isovector partner, the uq. We see from table 
4 that the width for the decay ijj(1650) ao7 is large, and from table 6 that it is a well-defined 
decay. So this decay can provide independent information on the existence and properties of the 
ao(1450). 



5.4 Hybrid meson radiative transitions 

Radiative transitions between hybrid states in the mass range under consideration are essentially 
ones in which spin-flip is required. As the qq pair is in a spin-singlet state for a 1 hybrid and 
in a spin-triplet state for ^, 1 ^ and 2 ^ hybrids then, analogously to (p^), for the transition 

(l^^)^f (J ^)h1 one has 

2J+1 4 Eb _2 .onN 

where Fh is a form factor similar to ( p^ ) , unknown but of the order of unity. If the mass ordering 
is indeed ^ < 1 ^ < 1 < 2 ^, the radiative decays of hybrid vectors to hybrid pseudoscalars 
and to exotic hybrids should be present but with widths less than 50 keV, the exact value depending 
on the phase-space. However if the mass ordering is the reverse, then the transition (0 
(1 )hJ is possible. With J = 1 in (|30|), then for the decay 7r//(1800) (jj/f(1500) this gives 
a rather healthy width of about 300 keV. Thus the decay 7r(1800) cj(1420) or cj(1650) is a 
potential discriminator among various assignments both for 7r(1800) and the two uj states. For 
example, if the 7r(1800) is a 3-^5*0 qq state, these radiative decays will be very strongly suppressed, 
with a width < 1 keV because of the orthogonality of the wave functions, unless the tj(1420) or 
w(1650) is 3'^ 5*1, which is highly unlikely. Equally, if the 7r(1800) is a 2^5*0 qq state, then these 
radiative decays will have a width of more than 1000 keV. Thus if a radiative width of several 
hundred keV is found, then the two states must be siblings, which should be most natural for 
hybrids. 

So we note a clear hierarchy. The radiative decay width is 0(1 MeV) for qq states in the same 
spatial state (25 to 2S or 35 to 35); it is 0(300 keV) for hybrid to hybrid; and it is 0(1 keV) 
for 35 to 25. Heuristically, the 25 to 25 is big because of phase space, the hybrid likewise but 
reduced by the gluon-quark spin coupling while the 35 to 25 is killed by the orthogonality of the 
wave functions. 

Radiative transitions between hybrids and qq mesons are the most uncertain. In constituent gluon 
models such transitions are highly suppressed as it is necessary to remove the gluon and to rearrange 
the colour degrees of freedom. In the flux tube model the mechanism is less explicit, and it has 
been argued that the flux tube is excited as readily as a quark, with no extra suppression for 
radiative transitions between hybrids and qq states. However, even with no such suppression, the 
radiative decay of hybrid vectors to '^Pj qq states is small as the transition is necessarily magnetic 
spin- flip and is also suppressed by phase space. This contrasts with the radiative decay of qq vector 
mesons to "^Pj qq mesons which proceeds in leading order by an electric transition. We shall report 
on this issue elsewhere. 
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6 Conclusions 



In table 8 we show the final states for the dominant decays listed in tables 4 and 5. Obviously the 
experimentally cleanest signals come from the tt+tt^ decay of /2(1270) in p''(1450) /2 (1270)7 
and the K~^K~ decay of /2(1525) in 0(1690) /2 (1525)7. The advantages are that the final-state 
mesons are comparatively narrow, their decays are two-body and there are no neutrals in the final 
state other than the photon. Additionally these are decays for which the calculation is well-defined 
and the answers reliable, as we discussed in relation to table 6. An experimental check on the 
validity of the quark-model approach, apart from the total rate, is provided by the decay angular 
distribution. From table 6 and (A. 6) we see that we expect the electric-dipole term to dominate 
and to give a nearly isotropic angular distribution dT/dcosO c>c 1 -I- cos^6'/13. Any significant 
deviation from this would imply an unexpectedly-strong contribution from the magnetic-dipole 
term. These two decays are unique identifiers of the p(1450) and the 0(1690) respectively. As 
discussed in the previous section, radiative decays of a hybrid pn or (pn to these final states is 
strongly suppressed and so there is no ambiguity. 



Decay 


Final state 


/9'^(1450) ^ /i (1285)7 


47r7, rj-KTT^ 


p"(1450) ^ /2(1270)7 


rnr'f 


w(1420) -> ai(1260)7 


7r+7r^7r'^7 


w(1420) a2(1260)7 




0(1680) ^ /i (1420)7 


KKtt^ 


0(1680) ^ /2(1525)7 


KKj 


0(1900) ^ /i (1420)7 





Table 8. Final states for the dominant radiative decays of tables 4 and 5. 

The decay 0(1900) /i(1420)7 discriminates between the /i(1420) and the 77(1440), as can be 
seen from table 5. The nearness of the masses and widths of these two states, and several common 
hadronic decay modes, have hitherto been sources of confusion. Although the magnitudes of 
the calculated radiative decays of the l^Di states are subject to some uncertainty, the width of 
0(1900) /i (1420)7 is necessarily much larger than that of 0(1900) — > r/(1440)7. In the latter 
case the overlap of the wave functions leads to a much stronger suppression of the decay than in 
the former case, in any model. 

The relative rates of the radiative decays of the p(1700) to /o(1370) and /o(1500), and of the 0(1900) 
to /o(1500) and /o(1710), change radically according to the particular model for qq-glueball mixing. 
The differences are sufficiently great, as can be seen in table 7, for the appropriate mixing scheme 
to be identified and the glueball mass determined. The physics that can be extracted from these 
decays is significant and merits every effort to overcome the experimental problem posed by the 
multiplicity of hadronic decay channels of the scalars. 

The observation of the decay Ci;(1650) ao(1450)7 is of high priority as it establishes consistency 
between the masses of the isoscalar and isovector nn states, with implications for the scalar glueball 
mass. Fortunately the ao(1450) has comparatively simple decay modes, such as KK. 

The larger partial widths should be measurable at the new high-intensity facilities being proposed. 
In some cases they may be measurable in the data from present experiments. We give two specific 
examples. The iorj decay of the a;(1650) has been observed in the E852 experiment |Q. If the 
w(1650) is the ID qq excitation of the w, then the "^Pq model gives the partial width for this decay 
as 13 MeV The partial width for the radiative decay Ci;(1650) ai (1260)7 is of the order 
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of 1 MeV, that is about 8% of the ujrj width. The E852 experiment has several thousand events 
in the iorj channel, so we may expect several hundred events in the 017 channel. Similarly both 
the /5{1450) and p(1700) are seen by the VES collaboration |Q in the prj channel with several 
thousand events. Both these states have strong radiative decays, the yo(1450) to /2(1270)7 and 
the /9(1700) to /i (1285)7. Assuming that the p(1450) and p(1700) are respectively the 25* and ID 
excitations of the p, then the "^Pq model gives the partial widths for the prj decays of the p(1450) 
and p(1700) as 23 MeV and 25 MeV respectively, so the radiative decays could again be present 
at the level of a few hundred events. 
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A Angular Distributions 



In the following equations the isospin factors / are given in (0) and (|), A in (|l2|) and /3 in ®. 
The form factors Fsi, Fs2, Fdi and Fd2 are defined in (|l6|), (20), ( p5| ) and ( p9| ) and the functions 
Ge and Gm in (|l|) and 
l^Si ^ l^So 



1 Eu , 9x9 

, n T.aV— At l-t- COS^ 61 

dcos6' 2 mA-m^- ' '^^ 



(A.l) 



23Si ^ l^So 



dV 

dcos t 



3i 



2AV 

12^ 



30 



dcos9 



(A.2) 



dT 

dcosd 



5/3" 



- 1 



4AV/ ^ 

32^2 



f^— -l) 



4AV V dVl 



3/32 V3/32/35 



9/32/32 



dcos9 



(A.3) 



23Si ^ l^Po 



dV 
dcosO 



-ap -{1 + cos^ ^)[^-^)[ 

9 TTiA mi \ 2/3^ / v 



El 

2/32. 



<2 _ P 



2/32 



16 V/32 



Gir]lF 



M S2 



(A.4) 



dT 
dcos 9 



= ap^^Ul - i cos2 9)Gl - :^(1 - cos2 9)GeGm 



24 V/3 



(^)"(1-cos2^)GL)/F|^ 



(A.5) 
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B Radiative Helicity Amplitudes 



In this appendix we show how the amphtudes used in the main body of this paper are related to 
helicity amplitudes that appear elsewhere in the literature. This also enables application of some of 
our results, which have been specified for e~^e~ annihilation, to be taken over to photoproduction. 

Following Gsl Eol the electromagnetic interaction may be written in the form 



Jg„ = '^ejfi(^ - 2is-px A- g^'^{k + k') ■ (B.l) 



where Cj is the quark charge in units of e, /i = gyja/2mq = O.lSGeV is the quark-scaled magnetic 
moment and A = t^J^i: /2uj e^^'^' . Ch oosing the z axis to be along the photon momentum, and 



e= —-^(1,1,0) for Jz = +1 then (BJ^) becomes 



Jem = 2J^iiY,ej{p{S^,+iSy)+g~\k^+iky)y'P' (B.2) 

V ^ j=l 



to be compared with (7.3) in |49 , whereby after summing over two quarks the matrix element 
becomes 



i\ -^i{e)ip{a+)RLo + g-'{L+)RLi) (B.3) 



where (e)^ = / of (m and 



Rli ^ (V'2i|e'^'-/'(-*)^|Voo) 

Rlo ^ (V^2o|e'^''/'lV'oo) (B.4) 

where the ipL.L^ are the bound-state wave functions for the for the qq state with orbital angular 
momentum L, |Q. 

Radiative widths are then 

^M:=pl Pl!ji±lu)^^iMS- ,B,5) 

and u! = \p\ for real photons; (e) = ^ for / = ^ / = 1 and (e) = ^ for / = 0(1) I = 0(1) and 
the Mx are helicity amplitudes with helicity A. 

B.l Ml transitions 

As an illustration and check on the normalisation we compute T^uj —^ ttj). 
The matrix element becomes ((cr+) = 1/V2) 

M+l = fly^RLO (B.6) 

for J] = +1 and (e) = 1/2. Hence 

3 

- |ap^4>oo(p)P (B.7) 
3 TO(j 4 
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to be compared with ( p^ with niq = 0.33, g = 1- In the harmonic oscillator basis 

2 

i?oob)^exp(-^) (B.8) 

which is the f3A = Pb limit of Fsi of (|l6|). Hence if /3 - 0.4GeV § ||, then 

T{uj 7r7) - QMIeV (B.9) 

In the text we show that this value for (3 does not fit weU with the detailed spectroscopy, and when 
realistic values are used, the actual width for w — > 7r7 is somewhat reduced. See figure 1. 

An analogous calculation for the radiative transition from the ijj{2'^Si) involves a radial wavefunc- 
tion and different magnitude for p, p* say. Thus 

^ .7) ^ V7 J 6 \w) p y^^j ^^-'^^ 

This also sets the scale for transitions between 7r(1800) and 6^(1420/1650) which potentially bear 
on the question of hybrid states. 

B.2 Application to excited states: El and M2 transitions 

Helicity amplitudes follow from the most general form of the single quark interaction with a trans- 
versely polarised photon in the algebraic form 

j;™ = {AL+ + Bcr+ + Ca,L+ + Da-L++)4ii^{e) (B.ll) 

In the non-relativistic limit C = D = and 

A{L+) = RLi/pg 

B{a+) = i?Lo/\/2 (B.12) 

We will be particularly interested in transitions between 25* and IP states. In the Gaussian wave 
function approach we have 



A{L+) = RLi/pg 

pgV^y 16/?V '^''P V 16/32 J - gpV3 

B{a+) = RL0/V2 

II3V2 -2 



= ^('-i^)'"ni6^)"^M''^^ ^^-''^ 

If we are only interested in terms up to 0{p'^) we can approximate the above as 



A{L, 



gpVs 
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The widths are then 

^9 



T{p{2S) fjj) = (2J+ E l^^^l' (B-15) 



A>0 

where A refers to the hehcity of the state and the matrix elements \M\\ are: 



As a specific application, 

T{pi2S) ^ /o7) 



/o: 


Mo = 


{A-B)/V3 


/i : 


Mo = 


A/V2 




Ml = 


{A-B)/V2 


/2 : 


Mo = 


iA + 2B)/V6 




Ml = 


{A + B)/V2 




M2 = 


A 


3V5 




\iA-B)/V3f 



(B.16) 



4 a 



pP^^lGE-fj^GM (B.17) 
27 mp V 4p"^ ' 



using /i/g = ^/a/2mq. The form in (|l^) reduces to this in the particular case where all /? values 
are the same and the isospin factor 7=1/4 

Note also the following translations: 

2^4^2 
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